Stochastic polarization formation in exciton-polariton Bose-Einstein condensates 
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We demonstrate theoretically the spontaneous formation of a stochastic polarization in exciton- 
polariton Bose-Einstein condensates in planar microcavities under pulsed excitation. Below the 
threshold pumping intensity (dependent on the polariton life-time) the average polarization degree 
is close to zero, whilst above threshold the condensate acquires a polarization described by a (pseu- 
dospin) vector with random orientation, in general. We establish the link between second order 
coherence of the polariton condensate and the distribution function of its polarization. We examine 
also the mechanisms of polarization dephasing and relaxation. 
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Introduction. — Bose-Einstein condensation (BEC) of 
exciton-polaritons in semiconductor microcavities has 
recently been demonstrated experimentally [l|, ]2j, y|- 
This important discovery opens the perspective for var- 
ious applications including room-temperature polariton 
lasers Q , microscopic optical parametric oscillators @ , 
and optical logic gates [6|. In this context, it is important 
to establish a reliable experimental criterion for the BEC 
of polaritons. From the point of view of the Landau the- 
ory of phase transitions jj > BEC requires the build up of 
an order parameter physically associated with the macro- 
scopic wave-function of the polariton condensate [8| . The 
phase of the order parameter is chosen by the system un- 
dergoing BEC spontaneously; the spontaneous symmetry 
breaking [9( is considered as the "smoking gun" for BEC. 
The order parameter build up is also accompanied by a 
decrease of the second order coherence parameter, <?2(0), 
indicating the change from a thermal distribution of po- 
laritons to a coherent distribution [101 ]. Unfortunately, 
due to the limited time resolution of the Hanbury-Brown 
and Twiss experimental set up, it is very hard to measure 



32(0) with a good accuracy [11 
In recent work 
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13l . Il4j it has been suggested that 



the build-up of the order parameter can be evidenced 
by polarization measurements. Effectively, the polariza- 
tion degree of light emitted by polariton condensates con- 
tains information on the amplitudes and relative phases 
of both components of the spinor wave function of the 
condensate. In this Letter we present the kinetic model 
of spontaneous formation of the polarization vector dur- 
ing the course of polariton BEC. Indeed, the polarization 
vector provides the same information on the condensate 
(complex) order parameter, as the second order coher- 
ence; this allows a straightforward method of evidencing 
the existence of spontaneous symmetry breaking in a po- 
lariton system under pulsed excitation. 

Existing theories for the dynamics of polariton BEC 
have been based either on the semiclassical Boltzmann 



equations describing the energy relaxation_of polaritons 
but neglecting the phase of the condensate 15], or on the 
Gross-Pitaevskii equation and its generalizations, assum- 
ing the existence of a coherent condensate from the very 
beginning [lfj, [lT], ll8(. The present model bridges the 
gap between these two approaches and allows one to de- 
scribe the formation of a coherent polariton condensate 
from an incoherent ensemble of polaritons. 

The formation of a polariton condensate is stochas- 
tic in its nature since polaritons entering the conden- 
sate have random phases and polarizations. Above the 
stimulation threshold when the average population of 
the lowest energy polariton state exceeds 1, the stochas- 
tic phase and polarization of the condensate start be- 
ing amplified and stabilize due to the stimulated scatter- 
ing of polaritons from an incoherent reservoir. Polariton 
interactions with acoustic phonons and between them- 
selves make the dynamics of the order parameter com- 
plex and non-trivial since they lead to dephasing and po- 
larization relaxation. Here we treat these effects within 
the Landau-Khalatnikov approach [19( generalized to de- 
scribe a two-component interacting Bose system. 

Formalism. — We consider BEC into one spin- 
degenerate ground-state level. The evolution of the 
order parameter if> a (t) is described by the Langevin type 
equation: 
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dt 
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Here a = ±1 denotes the two spin (circularly polarized) 
components of the order parameter [20(, W{t) is the in- 
come rate, T~ l is the polariton lifetime in the condensate, 
6 a (t) is the noise defined below. The third term in the 
righthand side of Eq. ([1]) describes the spin-dependent 
polariton-polariton interactions in the condensate 
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where a\ > defines the repulsion energy of two polari- 
tons with the same spin and «2 < defines the attraction 
energy of two polaritons with opposite spins. Finally, the 
last term in Eq. (Q]) describes the relaxation of the order 
parameter. The form of the functional $H rT [{ , (T ,-0*}] de- 
pends on the mechanism of relaxation. 

In the particular case of low condensate occupations, 
when one can neglect the polariton-polariton interac- 
tions and the order parameter relaxation, Eq. (fl]) reduces 
to the kinetic equation for the order parameter derived 
in 21| and originated from the quantum kinetic equa- 



tion describing the evolution of the condensate density 
matrix [22| . In this case the total intensity of the white 
noise, #o-(£), is given by the income rate of polaritons into 
the condensate, W(t). The correlator of the noise is 

{0 a {t)K'{t')) = \w{t)8^5{t-t'). (3) 

This noise is responsible for the phase and polariza- 
tion fluctuations in the ground state of the polariton 
system both below and above the condensation thresh- 
old. Above threshold, when the condensate occupation 
n(t) = IV'+il 2 + hA-i| 2 is n °t small, dephasing due to 
polariton self-interaction [23j becomes also important. 

In general, the income rate W(t) should be found by 
solving the semiclassical Boltzmann equation for the po- 
lariton relaxation into the condensate 11511. In what fol- 
lows, however, we adopt a simple model [171 ] considering 
all the polaritons that are not in condensate as a single 
incoherent reservoir. The reservoir occupation number 
N r (t) satisfies the kinetic equation: 



dN r 
~dT 



-T r N r -W(t)[n(t) + l] + P(t), 



(4) 



where P(t) is the incoherent pump rate and T" 1 is the 
life-time of polaritons in reservoir (and usually T r <C T c ). 
The exact dependence of the income rate on the reservoir 
occupation is defined by the relaxation mechanism. In 
the simplest case of polariton-phonon relaxation they are 
proportional to each other, W(t) — rN r (t). 

In the case of very short pulsed excitation, with pulse 
duration much less then r~ , the pump is reduced to 
the initial condition N r (0) — J P(t)dt for the reservoir 
concentration. Eq. (TJ| with P = is then solved simul- 
taneously with Eq. |(TJ) considering that the condensate is 
not initially populated, i.e., Vv(0) = 0. First, we present 
the results for the case when the polarization relaxation 
can be neglected (7 = 0), and then we describe the effects 
of relaxation of the order parameter. 

Formation of the order parameter. — The complex or- 
der parameter of the condensate ip a cannot be observed 
directly in a photoluminescence (PL) experiments. On 
the other hand, the polarization resolved PL gives access 
to the components of the condensate pseudospin (Stokes 
vector) [20J. These components are related to ip a as 

s x = (1/2)^;^.! + v+iV-ii), s y = (i/2)(V4 1 v_i - 
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FIG. 1: (a) The total polarization degree p of polariton con- 
densate as a function of pump intensity. Dots show the ran- 
dom values of p for different pulses, while the average (p) 
is shown by the solid line. The other solid line shows the 
second-order coherence (see text), (b) The distribution func- 
tion of the total polarization degree. Curves are labeled by 
the values of P/Pth ■ The time dependence of the average con- 
densate occupation number (n) is shown in the insert for the 
same values of P/P t h- 



V» + iV-i), S g = (1/2)(|t/; +1 | 2 - |^!| 2 ). Note that the 
condensate occupation, n, can be calculated from n 2 = 
4(S' 2 + S 2 + S 2 ). Averages and statistics of different ex- 
perimentally observable quantities are presented in this 
subsection. In what follows we will be interested mostly 
in time-integrated quantities which will be denoted by 
a bar, e.g., J ip a dt — tp a . The averaging over multiple 
pulses (i.e., over realizations of noise) will be denoted by 
angular brackets as in Eq. §3§ . In numerical calculations 
we used the parameters T c = 0.5 ps -1 , r r /r c = 0.01, 
ai = 1.8 MeV, a 2 /a x = -0.1, and r = lO^ps" 1 Q. 
The total polarization degree of the condensate, 



P = 



(^) 2 + (c-) 2 + (7-) : 



1/2 



(5) 



is shown in Fig. [TJ One can see that the solutions to 
Eqs. ([T]) and ([4]) exhibit a well pronounced threshold be- 
havior. The dynamical threshold 22j is defined by the 
balance of the polariton income and outcome rates for the 
condensate: W(0) = rN r (0) = r c , so that the threshold 
pump is defined by J Pthdt — T c /r. 



Below threshold, the average occupation of the con- 
densate is less than 1 and the order parameter fluctuates 
extensively during the reservoir lifetime T" 1 that defines 
the duration of the PL signal in this case. The pseudospin 
also fluctuates strongly both in amplitude and direction, 
so that the total polarization degree is close to zero, ft 
can be shown that (p) ~ (Tr/Tc) 1 / 2 for T r <C T c and 
P « Pth- 

Above threshold the condensate is formed during the 
formation time tf such that W(t) ^ T c for < t ^ tf, 
and disappears afterwards on the scale of T~ . This 
leads to a drastic increase in the condensate occupation 
and narrowing of the emission [see insert in Fig. [ljb)]. 
At the same time there is a strong increase in the 
total polarization degree of the condensate. Without 
polariton-polariton interactions the polarization degree 
would reach f . However, the polariton-polariton interac- 
tions [the third term in Eq. ([1])] result in a slow decrease 
of the average polarization degree for P > P t h ■ As a re- 
sult the dependence (p) on P exhibits a maximum placed 
close to the threshold. The build-up of the total polar- 
ization degree indicates the formation of a well-defined 
order parameter and well-defined pseudospin for each ex- 
citation pulse. It should be noted, however, that the di- 
rection of the pseudospin randomly changes from pulse 
to pulse. 

Formation of the order parameter at P > P t h can be 
observed also by measuring the second-order coherence in 
Hanbury-Brown and Twiss experiment [llj, [l2j, [25 1 . The 

second-order coherence parameter g( 2 )(0) = (n 2 )/(n) is 
also shown in Fig. [IJa). ft is seen that increase of the 
pump brings the polariton condensate from the thermal 
(chaotic) distribution with gW (0) = 3/2 for P < P th [II] 
to the coherent state with g( 2 > (0) ~ I for P 3> Pth, where 
the fluctuations are suppressed and the order parameter 
is well-defined for each pulse. Note that the polariton- 
polariton interactions that caused a decrease in the av- 
erage polarization degree for higher pump powers do not 
prevent <T 2 '(0) from approaching I. fn other words the 
polariton-polariton interactions do not result in loss of 
coherence; rather they have a dynamical effect on the 
pseudospin vector, which affects the average polarization 
degree as we shall discuss later. Clearly, both the second- 
order coherence and polarization measurements provide 
evidences of the condensate formation. However, the po- 
larization measurements give more details about the state 
of the condensate since the total polarization is sensitive 
to the effects of polariton-polariton interactions. 

Another important effect of polariton-polariton inter- 
action is the nonuniform distribution of the total polar- 
ization degree for P > P t h- Apart from strong fluc- 
tuations from to I present in this case the distribu- 
tion of polarization degree exhibits sharp and approxi- 
mately equidistant peaks that appear far above thresh- 
old as shown in Fig. QJb). The origin of these peaks can 



P/P th = 0.5 




P/P lh = 1.5 




PlP th =\.0 




P/P lh = 3.() 




FIG. 2: Showing the distribution function in normalized pseu- 
dospin space, s = 2S/n, for different values of the excitation 
pump. Far above the threshold the formation of the Larmor 
rings is seen. 



be understood if one examines the condensate polariton 
distribution function in the pseudospin space. This dis- 
tribution function is shown in Fig. [2] for four values of 
the pump. One can see that the in-plane component of 
pseudospin J S 2 + S 2 is strongly suppressed for several 

values of S z (see Fig. [3Jd) for P = 'iPth)- This happens 
as a result of spin-anisotropy of polariton-polariton in- 
teractions. Indeed, the interaction Hamiltonian |(5J) can 
be re-written as H int — (l/4)(ai + a^rv 2 + {a,\ — a.-i)S 2 zl 
and the third term in ([I]) produces the self-induced Lar- 
mor precession [27J of pseudospin around z-axis. The 
frequency of this precession is proportional to S z . For a 
fixed n the in-plane component of pseudospin is averaged 
to zero for the values of S z that give full 2ir rotation, fn 
reality the condensate occupation is not fixed and fluc- 
tuations also exist, so that the in-plane pseudospin never 
averages to zero exactly. Nevertheless this effect gives rise 
to a sequence of "Larmor" rings seen at strong pumping. 

Relaxation of the order parameter. — The results ob- 
tained above are valid for the case of short condensate 
lifetime T" 1 compared to the typical relaxation times of 
the order parameter, fn this subsection we show that fast 
relaxation modifies qualitatively the polarization proper- 
ties of the polariton condensate. Relaxation not only 
smears out the Larmor rings discussed above but also 
favors a linear polarization of the condensate. 

Below we consider the simplest form of the relaxation 
term, namely, model A also referred to as the Landau- 
Khalatnikov or Onsager model [l9(. In this model, the 
relaxation term in Eq. {1} is written as ^R CT = SHint/Sip^. 
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FIG. 3: Showing the effect of relaxation of the order param- 
eter on distribution function in normalized pseudospin space, 
s = 2S/n. Contrary to the case without relaxation (Fig. [2]| a 
linear polarization is formed for strong pump above threshold. 
The relaxation parameter 7 = 1 meV -1 ps _1 . 



This term favors relaxation of the order parameter to 
minimize the polariton-polariton repulsion energy. This 
repulsion energy has its minimum if half of the polaritons 
in the condensate have their spins oriented up, whilst 
the other half have their spins down, corresponding to a 
suppression of the z component of the pseudospin vector 
and the formation of a linearly polarized condensate. It 
should be noted that the relaxation term we consider does 
not conserve the number of polaritons in the condensate. 
In this way, the condensate depletion and leakage of po- 
laritons from the condensate are implied by the present 
model. 

The distribution functions of polaritons in the normal- 
ized pseudospin space for different pump intensities are 
shown in Fig. [3l One can see that the circular polar- 
ization is rapidly lost by the condensate as the pumping 
intensity increases. The characteristic oscillations of the 
average linear polarization of the condensate as a func- 
tion of its circular polarization degree seen in Fig. [2 van- 
ish once the relaxation of the order parameter is switched 
on. The distribution of polariton condensates in the 
pseudospin space changes from a sphere to torus and fur- 
ther approaches a flat ring shape as the pumping power 
increases. The condensate develops a strong spontaneous 
linear polarization at strong pumping. The build-up of 
the linear polarization as a result of the polariton BEC 
has been experimentally observed [l|, |2j, while in these 
works the orientation of linear polarization plane has 
been pinned to one of the crystal axes of the structure due 



to some intrinsic anisotropy of the samples [l| . Further 
experimental research is needed in order to evidence the 
stochastic vector polarization of a polariton condensate. 

Conclusion. — We have described theoretically the 
spontaneous build-up of the vector order parameter in 
the course of polariton BEC. The build-up of the order 
parameter manifests itself in formation of the stochastic 
vector polarization of the condensate. This polarization 
is correlated with the second order coherence of the con- 
densate. Its statistical distribution depends on the mech- 
anism of polarization relaxation in the polariton system. 
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